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1 Abstract 


"This is the version of the article before peer review or editing, as submit- 
ted by an author to the Astrophysical Journal. IOP Publishing Ltd is not 
responsible for any errors or omissions in this version of the manuscript or 
any ion derived from it. The Version of Record is available online at 
[https://doi.org/10.3847/1538-43577/abd4 ec] ". 

As part of the NASA Starlight collaboration, we look at the implicatioi 
radiation effects from impacts with the interstellar medium (ISM) on a di 
energy driven relativistic spacecraft. The spacecraft experiences a stream of 
MeV /nucleon impacts along the forward edge primarily from hydrogen and he- 
lium nuclei. The accumulation of implanted slowly diffusing gas atoms in solids 
le proc s of bubble formation, blistering, 
and exfoliation. This results in macroscopic changes to material properties and, 
in the cases of blistering and exfoliation, material erosion via blister rupture and 
delamination. Re stic hydrogen and helium at constant velocity will stop in 
the material at a similar depth, as predicted by Bethe-Bloch stopping and subse- 
quent simulations of the implantation distribution, leading to a mixed hydrogen 
and helium system similar to that observed in fusion plasma-facing components 
(PFC’s). However, the difference in location of near-surface gas atoms with 


drives damage through the meso- 


respect to the direction of exposure means that previously developed empirical 
models of blistering cannot be used to predict bubble formation or blistering 
onset. In this work, we present a model of the local gas concentration threshold 
for material blistering from exposure to the ISM at relativistic speeds. Expected 
effects on the spacecraft and mitigation strategic " 
considerations apply to the Breakthrough Stars 


2 Introduction 


2.1 Motivation 


1f exploration beyond the outer limits of the solar system is to be conducted 
within reasonable fractions of a human lifetime, a radical shift in spacecraft 
propulsion technology is necessary. The radical shift proposed by Lubin [Lubin(2016)] 
and the NASA Starlight project consists of the directed energy propulsion of 
ultra-light spacecraft and a low absorption reflective sail structure. The baseline 
system propels a low mass spacecraft via photon momentum exchange using a 
1-10 km aperture, 100 GW class laser phased array to relativistic speeds. Such 

a system can be used with a wide variety of missions including the acceleration 

of a 1 g spacecraft to ~c/4 or a 10 kg spacecraft to ~c/40 or greater. This 
would enable the first robotic interstellar journey to the nearest star system, a 
Centauri, with a flight time of less than 20 years. For this work, we use the term 
relativistic spacecraft to refer to this design and mission velocities greater than 
0.1c. A spacecraft travelling at these speeds through the ISM would experien 1 

the ISM in its rest frame as a nearly mono-energetic beam of particles 
ing primarily of protons, electrons, and alpha particles. Thus, in designing a 
spacecraft capable of surviving interstellar journeys, it is crucial to character- 
ize, understand, and mitigate the damage that is caused by incident particles 
to ensure spacecraft survival through the ISM. 


2.2 ISM Composition and Homogeneity 


The ISM is composed of roughly 99% gaseous matter and 1% granular dust 
by mass. Of the gaseous matter, approximately 70% is hydrogen and 28% is 
helium by mass, while the remaining 2% is made up of heavier elements including 
carbon, oxygen, and iron [Klessen & Glover(2016)] [Draine(2011)]. Interstellar 
dust is composed of still heavier constituents including hydrocarbons, silicates, 
and ices in the form of crystalline and amorphous grains of size less than ~1 jim. 
[Draine(2011). Hydrogen in the ISM is observed in both neutral and ionized 
states, where the neutral state constitutes atomic and molecular hydrogen. The 
gaseous component of the ISM is commonly divided into five dominant phases 
which, since the phases are generally not in thermal equilibrium, drive large 
scale dynamics that produce structures like interstellar gas clouds of varying 
temperatures and often promote smaller scale gravitational phenomena such as 
star formation. 


On large scales (~10 to 100 ly), the ISM is inhomogeneous and can vary 
greatly in density. While the majority of the ISM by volume is composed of 
ionized hydrogen with a density of approximately 1 cm-?, most of the neutral 
hydrogen is collected in large, cool clouds with densities between 0.1 and 50 
cm? for neutral atomic hydrogen clouds and often greater than 1000 cm~? 
for clouds cool enough to allow molecular hydrogen to exist in their cores 
Klessen & Glover(2016)] [Draine(2011)]. These neutral gaseous phases are com- 
monly referred to as cold and warm neutral media (CNM and WNM, respec- 
tively) depending on their temperatures, while the ionized phas 
into the warm and hot ionized media (WIM and HIM, respectively). 

For relativistic spacecraft for which a Cen or any other local star is the 
target, the ISM can be approximated by measurements of the local interstel- 
lar cloud (LIC) due to the relatively short distance from the Sun to a Cen- 
tauri. Common measurements of LIC densities, such as those performed by 
Gloeckler and Geiss, report combined H and He densities as low as ~0.26 cm? 
Gloeckler & Geiss(2004)]. Therefore, for our analyses we conservatively ap- 
proximate the local ISM as a homogeneous gas of atomic hydrogen with a density 
of 1 cm? and a temperature of approximately 8,000 K. Thus, to a relevant first 
approximation, the local ISM falls within the WNM phase. However, distinct 
structures produced by the interaction of the ISM and stellar winds complicate 
the environment since they can create phases of matter that deviate greatly 
from the above approximation. 


2.3 Stellar Winds 


In addition to the various phases of the ISM that the spacecraft may experience, 
it may also be subject to the potentially damaging effects of stellar winds and 
the structures they produce as it approaches stellar targets. Specifically, in the 
simple stellar wind model consisting of a radial and isotropic outflow of material 
from the star, one would expect that, in the region where the stellar wind inter- 
acts strongly with the local ISM, there would exist a potentially different phase 
of material created by the interacting materials. In a NASA study by Linsky 
and Wood [Liusky & Wood(1996)], it is described how the a Centauri sy 
though consisting of a total of three stars, likely produces a solar wind and astro- 
sphere similar to the Sun’s due to the proximity of the in the system, the 
main sequence nature of the two dominant stars, and the similarity in relative 
velocity to the local ISM of the Sun and the a. Centauri system. It was found, 
and later verified by a number of groups [Zank et al.(2013)| [Wood et al.(2004)], 
that in the region surrounding the a Centauri astropause, and thus likely also 
around the heliopause of the Sun and the astropause of most stars in the local 
cloud, there exists a heated region of ISM material compressed by the outflow- 
ing stellar wind to a temperature ~30,000 K and density ~0.3 cm-?, commonly 
referred to as the hydrogen wall. It is interesting to note that at 0.25c, the sp 
craft will likely pass through the hydrogen wall very quickly. Thus, we may be 
able to test this regime experimentally since the exposure time at hydrogen wall 
and temperatures may be small. 


em, 


Figure 1: a) Prototype 100 mm diameter silicon wafer scale spacecraft with a 
plasma etched honeycomb pattern to reduce the mass to 1 g. b) Detail of an 
etched honeycomb vertex. Fabricated at CNSI/UC Santa Barbara. 


2.4 Relativistic Spacecraft Design 


The baseline design for small relativistic wafer scale spacecraft consists of a 


single thin semiconductor or hybrid disk of radius between 25 and 100 mm 
depending on the mission with a nominal diameter of 100 mm [Lubin(2016)). 
The spacecraft flies edge on to reduce the number of ISM particle and dust 
impacts. The wafer houses the entire spacecraft bus, sensor suite, power system, 
and communication system on its surface, and would be manufactured from 
silicon or titanium using plasma etching technologies. It is to be extremely thin, 
~100 um, with a thin-walled honeycomb structure on the other side to reduce 
mass and enhance its structural integrity, as shown in Figure[I] The spacecraft 
wafer is then embedded within or attached to a light-sail structure made of an 
ultra-thin (~0.1 to 1 jm) dielectric, reflective and low absorption material that 
reflects the incident laser light, thus generating thrust. In other variants, the 
spacecraft could be a long thin “needle-like” structure with small forward-edge 
cross section. The lowest mass missions have a wafer and reflector mass of order 


1 g, though the same launch system is capable of larger missions with higher 
scraft being slower with speed scaling as m~!/* [Lubin(2016)]. 

In some proposed mission scenarios, the spacecraft design calls for a mec 
anism by which the reflector can be ejected from the wafer after the spacecraft 
has been accelerated to its cruise speed. In others, the reflector remains at- 


tached to the spacecraft for the duration of its journey or is removed passively 
by vaporization due to collisions with interstellar gas and dust. 


Due to the extremely large distances to interstellar mission targets, the 
spacecraft communicates with Earth using a dedicated laser communications 
stem to send data back to Earth. The spacecraft quickly enters a regime 
where the command delay time becomes exceedingly long, requiring largely 
The spacecraft is 
powered by the heat given off by a small pellet of radioactive material such as 


S; 


autonomous operation with infrequent uplink command; 


plutonium-238, though the ISM impacts could conceivably be used as a source 


of power. In most mission scenarios, the spacecraft enters a low power hiberna- 
tion mode during transit to its destination. The details of this are discussed in 


our other papers [Lubin(2016)]. 
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Figure 2: Schematic of the various particle interactions in the plane of the 
spacecraft upon impact of ISM protons and heavier species on its leading edge. 
With exaggerated scales, the shield is shown in grey and the spacecraft wafer 
in blue. As shown, while the protons are stopped by the shield at the H/He 
Bragg peak, secondary and tertiary particles as well as heavier species have 
the capability of impinging far deeper into the spacecraft. 


3 Cumulative Effects of ISM Impacts 


3.1 Instantaneous Ion-Matter Interactions 


ISM particle impacts will have a significant effect on relativistic spacecraft. 
Figure []shows a summary of particle-spacecraft interactions, including hydro- 
gen/helium implantation, heavy species impacts, and secondary particle genera- 
tion. Since ISM particles incident upon a relativistic spacecraft will immediately 
lose their electrons on impact, we can treat the ISM in the frame of a relativis- 
tic spacecraft as a fully ionized, constant-velocity beam, composed primarily of 
hydrogen and helium. Ion-matter interactions span many orders of magnitude 
of energy-, time-, and length-scales and include the processes of sputtering, dis- 
placement of target atoms and damage, thermal spikes, and implantation. Con- 
temporary understanding of these phenomena comes from a combination of the- 
oretical, empirical, and computer models. Careful application of these models is 
necessary to understand the cumulative effects of ion bombardment. For exam- 
ple, a straightforward calculation of the analytic Bohdansky [Bohdansky(1984)| 
or semi-empirical Yamamura [Matsunami et al.(1984)] formulas for the sputter- 
ing yield of relativistic ISM impacts will predict a negligible amount of material 
erosion during an interstellar journey, as shown in Figure[}] (a). 

Total one-dimensional erosion depths can be calculated based on the fluence 
and target material parameters with the expression 
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Figure 3: (a) Sputtering yields from the Yamamura and Bohdansky formulas 
for hydrogen and helium incident on copper and silicon. (b) Total erosion from 
the Yamamura formula for sputtering yield, expected ISM fluence, atomic 
mass of the shield material, and mass density of the shield material, as shown 
in Equation [I] 


where Az is the depth of material eroded, 9; is the fluence of species 1, Yi2(E1) 
is the sputtering yield for species 1 incident on species 2, Mo is the atomic mass 
of species 2, and p» is the mass density of species 2. Using this approxima- 
tion, erosion of silicon and copper by hydrogen and helium for a journey to 
a Cen at speeds between 0.1 and 0.5c is less than a single atomic layer, as 
shown in Figure [3] (b). Erosion by physical sputtering is so low for two rea- 
sons: first, energy transfer during collisions between energetic light ions and 
heavy atoms is inefficient; second, most of the ion-atom energy transfer occurs 
deep below the surface, where the nuclear stopping cross section becomes sig- 
nificant compared to the electronic stopping cross section. For atoms to be 
sputtered, momentum must be transmitted from this location to a free surface. 
A previous work has considered the damage caused by the 
cumulative effect of instantancous ion-matter interactions, including sputtering 
and track-formation. However, damage from the accumulation of gas atoms is 
an unexplored, potentially threatening phenomenon for relativistic spacecraft. 
Gas atom accumulation negatively affects material properties, causes material 
swelling, induces surface morphology changes through blistering, and can drive 
erosion through exfoliation at rates exceeding physical sputtering. 


3.2 Bubble Formation and Blistering 


Bubble formation, blistering, and exfoliation are deleterious phenomena caused 
by the accumulation of relatively insoluble gas atoms in solids. Figure [7] shows 
an illustration of the bubble formation and blistering process. Insoluble gas 
atoms implanted in a solid, such as helium, can aggregate and form overpressur- 
ized bubbles [Donnelly (1985)]. If the local gas concentration exceeds a critical 


Figure 4: Diagram showing damage to materials in the intermediate fluence 
range, including material defects, blistering, and exfoliation, shown here for 
illustrative purposes. Reprinted by permission from Springer-Verlag: 

Scherzer(1983)}, 1983. 


concentration, near-surface bubbles can deform the surface and burst, in pro- 
s called blistering and exfoliation. When blistering does not occur, such 
as in the case of a subcritical gas concentration, the implanted gas atoms neg- 


c 


atively affect material properties and can lead to swelling and crack nucleation 
Condon & Schober(1993)]. Figure [I] shows examples of surface damage caused 

by implanted gas atoms, including blistering and exfoliation. 

Blistering is a phenomenon typically encountered when materials are ex- 

posed to high-flux, high-fluence, intermediate energy (keV range), insoluble gas 

ions, such as in ion beam experiments and nuclear fusion [Balden et al.(2013) 

Behrisch(2010)] [Bauer(1978)]. Blistering by hydrogen and helium specifically 

has been extensively studied in these contexts [Das & Kaminsky(1975)] [Kaletta(1980)] 

St-Jacques et al.(1978)] [Evans(1976)] [Erents & McCracken(1973)] [Guseva & Martynenko(1981)]. 
Blistering induced by low energy protons has also been observed in extreme ul- 

traviolet multilayer optics used in photo-lithographic systems 

and space instrumentation [Pelizzo et al.(2011)|, and in metal films tested for 

space environmental conditions [Sznajder et al.(2018)]. n example 

of blistering occurring in a thin gold layer is shown. The sample has been irra- 

diated with 4 keV low energy helium ions with a total fluence of 4x 10?! m-?; at 

lower fluences, no clear evidence of bubbles was found. At higher energies (MeV 


range), blistering has been reported to occur at fluences an order of magnitude 
below the critical dose found from analogous experiments at intermediate en- 
ergies (keV range). [Gavish Segev et al.(2017)]. In these situations, blistering 
occurs on the surface exposed to the ion flux, because the ions implant close to 
the surface. 

In the frame of a relativistic spacecraft, the ISM appears as a low-flux, high- 
fluence, MeV-energy beam that spans the entire front-facing cross section of the 


Figure 5: (a) Gold sample prior irradiation. (b) Gold sample after irradiation 
with 4 keV He ions with a total fluence of 4 x 10'7 cm~?. Blistering is observed 
close to the surface at ions penetration depth. Reprinted with permission from 
- Copyright (2018) American Chemical Society. 


spacecraft. During relativistic travel, the implantation depth distributions of 
ISM gas atoms will be strongly peaked many atomic layers below the surface. 
Additionally, due to the scaling of Bethe-Bloch stopping, at a constant relati 
tic velocity, hydrogen and helium will stop at approximately the same depth. 
From the prefactor of the Bethe formula, So, shown in Equation]with Zi being 
the atomic number of the incident ion, Z2 being the atomic number of the target 
material, m, being the electron mass, and v being the incident velocity, helium 
will experience a stopping power approximately 4 times larger than hydrogen. 
However, since helium is approximately 4 times more massive than hydrogen, 
helium at constant velocity will have approximately 4 times more kinetic energy, 
resulting in similar ranges in a given material. This will result in a thin layer of 
material spanning the entire exposed portion of the spacecraft, at the depth of 
the implantation distribution, where the local concentration of mixed gas atoms 
may form bubbles and may exceed the critical concentration necessary to form 
blisters. The uniformity of the implantation across the entire front-facing cros 
section suggests near-surface bubbles will form almost immediately on surfaces 
perpendicular to the ion flux. 
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In Figure [6] (a) 
beam experiment is 
of bubble formation is approximately the range of the implantation distribution. 
In order to cause blistering or modify surface morphology, bubbles must travel 
to this surface to have an effect. Figure [6] (b) shows bubble formation caused 


an illustration of bubble formation during a high energy ion 
own; the distance to the nearest surface from the location 


ISM in spacecraft fame 
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Figure 6: (a) Bubble formation due to high energy ion beam exposure with 
beam width significantly narrower than target size. Rs, the distance from 
bubbles to the nearest surface, is close to the range of the ion implantation 
distribution, approximated above by a normal distribution. (b) Interstellar 
spacecraft configuration traveling edge on. Rs is effectively zero because ISM 
exposure is spacecraft-spanning. Any swelling, blistering, or exfoliation will 
likely occur first on surfaces with normals perpendicular to the incident flux. 


by ISM gas accumulation. In this case, the smallest distance from bubble for- 
mation locations to a surface is effectively zero; bubbles do not need to travel 
any distance to begin having a deleterious effect on the nearest surface. For this 
reason, swelling, blistering, and exfoliation may happen significantly earlier for 
the case of wide exposure. Previous work on blistering has produced an empiri- 
cal formula for predicting the dependence of critical 
dose on the incident ion energy; since this formula assumes blister formation 
at the surface directly exposed to ion flux, it is not applicable here. Instead, a 
more fundamental model of blistering onset must be used. 

A simple theoretical model of blistering onset was described by Martynenko, 
and a reduced version is reproduced here [Martynenko(1977)]. In this model, a 
local pressurization is determined from the concentration of dissolved gas atoms, 
cg, and the dissolution energy, H, of the gas species in the target material: 


p= Hcg (3) 


A critical pressure for blistering onset is chosen as the pressure which exceeds 
the appropriate yield strength of the material, øy, resisting the overpressur- 
ization in implanted gas bubbles. A critical concentration, ce, is the local gas 
concentration at which this pressure is exceeded. Local gas concentration is 
approximated by: 


p 
7 AR (4) 
where © is the fluence and AR is the width of the distribution. The width of 
the distribution is determined from the width of the implantation distribution 
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Figure 7: Illustration of the process of bubble formation and blistering, 
neglecting gas atom interactions and material defects. a) At very low fluence 
(< 1075 m-?), individual gas atoms are implanted in the material. They may 

diffuse freely and aggregate, or be trapped by defects. b) At low fluence 

(71075 m-7?), nano-scale bubbles will form in the material, but no surface 
effects will be visible. c) At intermediate fluence (1015-10?! m~?) the local gas 
concentration will put stress on the material, and if at the surface, will begin 
to make morphological changes. d) At high fluence (> 102! m~?) surface 
blisters will begin to burst, leading to exfoliation, flaking, and ultimately, 
material erosion. 


and the diffusion length: 


AR = AR? + Dt (5) 


Where AR? is the straggle or standard deviation of the implantation distri- 
bution, D is the diffusion coefficient, and t is the irradiation time. In the 
case of gas implantation via ion irradiation, it is often assumed that irra- 
diation damage produces enough defects to serve as trapping sites that the 
diffusion coefficient is negligible [Martynenko(1977)]| [Scherzer(1983)]. Lattice 
defects caused by ion irradiation provide nucleation for fixed helium bubbles 


Kornelsen(1972)]. We can assume additional trapping due to the interaction of 


hydrogen with helium in the material, since helium bubbles serve as trapping 


sites for hydrogen atoms through synergistic effects [Hayward & Dco(2012)], 


leading to a lower critical dose in the case of mixed hydrogen and helium ex- 
posure [Guseva & Martynenko(1981)]. At temperatures well below the melting 
point, such as that expected for a relativistic interstellar spacecraft, high-flux hy- 
drogen irradiation damage effects have been reproduced in low-flux experiments 
[Gao et al.(2019)]. However, for an interstellar probe, the flux is sufficiently low 
that diffusion effects may play a non-negligible role. Classical diffusion of indi- 
vidual gas atoms could lead to lower local gas atom concentrations as the implan- 
tation distributions widen, or increased gas atom concentrations around trap- 
ping sites such as grain boundaries. Helium bubble nucleation and migration 
could lead to increased gas atom concentrations and damage near surfaces, grain 
boundaries, and other defects [Nakamura et al.(1977)] [Lane & Goodhew(1983)] 
Goodhew(1983). No single theoretical framework exists to summarily treat 
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Table 1: Range and straggle of hydrogen and helium at 0.2c 


Material | C (graphite) | | Cu Si 

Ru(0.2c) T9mm |0.72mm | 2.1mm 

Rue (0.2) T9mm |0.7Imm | 21mm 
AR? (0.2c) 30jan Ibam | 354m 
AR? s. (0.2c) 13um T2um | lym 


these effects. However, simple models such as that presented here immediately 
offer compelling mitigation strategies. To perform this analysis, we find im- 
plantation profiles of ISM gas atoms at relativistic speeds using a BCA code, 
calculate critical concentrations for blistering onset for hydrogen and helium 
individually assuming a wor ario of negligible diffusion, and show the 
effect of non-negligible diffusion on local gas concentrations. 


3.3 Ion Implantation Distributions 


SRIM is a free-use but closed-source Monte Carlo, Binary Collision Approxi- 
mation (BCA) code used to model ion-solid interactions in continuous devel- 
opment since 1985 [Ziegler et aL.(2010)] [Ziegler(2004)]. A detailed description 
of the BCA model can be found in Eckstein [Eckstein(1991)| and Robinson 
[Robinson(1994)|. For nuclear interactions, SRIM uses a universal interaction 


potential, the Ziegler-Biersack-Littmark potential, and the MAGIC algorithm to 
calculate the scattering angle of each binary collision [Biersack & Haggmark(1980)]. 
To model electronic interactions, SRIM includes an electronic stopping formula 
that includes Lindhard-Scharff [Lindhard & Scharf(1953)] electronic stopping at 
low energies (below approximately 25 keV/amu), Bethe-Bloch [Ziegler(1999) 
stopping with corrections at high energies (above approximately 1 Mev/amu), 
and Andersen-Ziegler [Andersen & Ziegler(1977)] stopping at intermediate en- 
ergies between the ranges of validity of Lindhard-Scharff and Bethe-Bloch. Elec- 
tronic stopping is evaluated along each path between binary collisions with tar- 
get atoms. SRIM produces detailed information about stopped ions, from which 
we find implantation distributions. 

Although SRIM has recently come under scrutiny for inaccuracies for low to 
intermediate energy sputtering (Shulga(2019)] [Shulga(2018)] [Wittmaack(2004)] 
and sputtered atom angular distributions [Hofsáss et al.(2014)], these issues do 
not affect the implantation distributions at high energy used in our ana 
From these ion implantation distributions we calculate the straggle, o: 
dard deviation of the implantation distribution, and use this to estimate the 
local gas concentration in materials. For this calculation, ISM partic 
sumed to be mono-energetic for the duration of the mission, and the acceleration 
phase is not considered. 

Figure | (a) shows the implantation depth and straggle of hydrogen and he- 


an- 


s are as- 
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Range and Straggle in Cu 


(a) 


Depth Datrbuors in Cu 


Figure 8: (a) Implantation depths for H and He in copper from SRIM 
simulations. The regions bounded by solid lines for hydrogen and dashed lines 
for helium are centered on the mean implantation depth and span above and 
below the mean by the straggle, +AR. (b) Implantation distributions, 
normalized and scaled to the appropriate ISM species number density, for H 
and He in copper from SRIM simulations at 0.1, 0.2, 0.3, 0.4, and 0.5c. 
Hydrogen distributions are unfilled and bounded by solid lines, and helium 
distributions are filled and bounded by dashed lines. Heights are normalized to 
1 for all distributions. 
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lium in copper, an example shield material, from 0.1 to 0.5c. Copper was chosen 
as a representative metallic shield material in which hydrogen is relatively solu- 
ble, potentially mitigating hydrogen blister formation [Magnusson & Frisk(2017)]. 
For velocities in excess of approximately 0.07c, the implantation distribution of 
helium out to one standard deviation from the mean fits entirely within one 
standard deviation of the hydrogen distribution. Inset is a zoomed in figure 
showing this overlap at velocities relevant to a relativistic, interstellar mission 
(0.25-0.3c). Any long duration, high velocity mission will be subject to a layer 
of mixed hydrogen and helium implanted many atomic layers deep below the 
surface. Figure [B] (b) shows the full implantation distributions of hydrogen 
and helium in copper at 0.1, 0.2, 0.3, 0.4, and 0.5c from SRIM simulations, 
normalized to equal height. These results show that at speeds relevant to an 
interstellar mission, we expect a layer of mixed hydrogen and helium implanted 
in the spacecraft. Table[I]shows the results of SRIM simulations for a number of 
other example materials at 0.2c. In each material, hydrogen and helium stop at 
me distance. As the velocity increases, the straggle increases, due to the 
sed number of small-angle atomic collisions incident particles are subject 
to as they travel deeper in the material. 


3.4 Determination of Blistering Onset 


A relativistic spacecraft on a mission to a Cen will likely be travelling between 
10 and 30% the speed of light. In the frame of the spacecraft, hydrogen from 
the interstellar medium will be impacting with an energy of 4.7-45.3 MeV, and 
helium will be impacting with an energy of 18.8-180.0 MeV. Assuming an av- 
erage ISM number density of 1.0 cm-?, with the hydrogen contributing 92% 
and the helium contributing 8% by number, the resulting flux for each species 
at 0.2c will be 2.8 x 10!? s-!m-? and 2.4 x 107? s-'m-? for hydrogen and 
helium respectively. Over the course of a journey of 4.37 lyr, the approximate 
distance to a Cen, the fluence for hydrogen and helium will be 3.8 x 10? m-? 
and 3.3 x 10?! m~? respectively, regardless of travel speed. 

Using the model for determining blistering onset presented above, Figure [f] 
shows calculated critical concentrations of implanted gas atoms for a range of 
dissolution energies. Hydrogen, which dissolves readily in copper, has a di: 
lution energy near 0.5eV [Magnusson & Frisk(2017)]. For gas-material comb 
nations with lower solubility, the dissolution energy will be significantly higher. 
Dissolution energies of hydrogen and helium are not widely available for many 
materials. Determining dissolution energies requires carefully designed experi- 
ments or computational modeling. Additionally, this calculation assumes a con- 
stant ISM number density of 1.0 cm-?; if the actual number density is larger, 
the fluence and thus implanted gas concentrations will increase linearly with 
increasing number density. From Figure [J] it is apparent that the local gas 
concentration of hydrogen and helium at 0.1c will be lower than at 0.2c. At 
higher relativistic veloc incident ions will have travelled through signifi- 
cantly more material and will have experienced many more small-angle nuclear 
collisions, leading to increased straggle. 
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Figure 9: Critical implanted gas concentrations for blistering onset calculated 
for hydrogen and helium assuming ISM particle density on the order of 1.0 
cm-? for dissolution energies of 0.1, 1.0, and 10.0 eV in a material with yield 
strength oy. Horizontal lines marking estimated gas concentration in copper 
for hydrogen and helium at 0.1 and 0.2c are shown for reference. Vertical lines 
marking estimations of the yield strength against blistering in copper are 
shown. 250 MPa is the normal yield strength of copper. 13.3 GPa is 10% of 
the elastic modulus of copper, per [Martynenko(1977). 


At an ISM density of 1.0 cm-?, neglecting diffusion and synergistic hydrogen- 
helium effects, and using the straggle of hydrogen and helium in copper, this 
model suggests that materials with a yield strength well in excess of 10 GPa will 
not blister. Determining the appropriate yield strength presents some difficulty. 
Using standard yield strength significantly underestimates the critical fluence 
for many materi one egy is to estimate the yield strength as 1096 of the 
elastic modulus Values for these quantities are reported for 
copper in Figur t 250 MPa and 13.3 GPa, representing the 
traditional yield strength h of the elastic modulus respectively. Since we 
expect near-surface bubbles on surfaces perpendicular to ion exposure, choosing 
the correct yield strength will be critical for choosing target materials. This 
value may be obtainable from molecular dynamics simulations of near-surface 
bubbles. Note that this figure shows the threshold for blistering, and that bubble 
formation, swelling, and changes to material properties are expected at fluences 
many orders of magnitude lower than the blistering threshold, and these effects 
of gas accumulation will need to be mitigated. 


3.5 Role of Diffusion in Cumulative Damage 


When our analysis extends to the case of significant classical gas species dif- 
fusion, we can determine for what diffusion coefficients diffusion will have a 
measurable effect on the local gas concentration. Figure[I0|shows the estimated 
local gas concentration in a copper shield for a spacecraft traveling at 0.2c. Dif- 
fusion is only significant when the diffusion length becomes comparable to the 
straggle of the ion distributions. In copper, classical diffusion does not become 
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Gas Concentrations in Capper with Difin 


Figure 10: Gas concentration in copper after a journey to a Cen at 0.2c. 
Diffusion is only significant when the diffusion length becomes comparable to 
the width of the implantation distribution; in copper, this occurs when the 
diffusion coefficient is near 10-17 m?s-!. This is on the order of classical 

diffusion coeffi 


nts in metals. 


ra 


significant until the diffusion coefficient exceeds 10717 

This analysis assumes an isotropic diffusion coefficient for gas species, but 
implanted gas atoms can have significantly different transport behavior. For 
example, helium atoms implanted in many materials aggregate and form bub- 
bles. Helium and hydrogen bubbles are known to have a synergistic effect in 
driving deleterious effects in solids [Hayward & Deo(2012). This interaction 
can lead to increased local gas concentration, amplifying the effect of bubbles 
on the material. 


Tn a radiation damaged material, insoluble atoms may interact strongly with 
material defects such as vacancies, significantly limiting diffusion. For example, 
in body-centered cubic (BCC) metals such as tungsten, a candidate material 
for fusion reactor PFCs, the hydrogen diffusion activation energy 
on the order of 0.1 eV, while the binding energy of a single va 


, cas 
been 
a signifi- 
atoms in 


pable of trapping as many as 10 atoms of hydrogen simultaneous 
calculated to be as high as 1.2 eV [Johnson & Carter(2010)]. This 
cant barrier to diffusion. Understanding the behavior of implanted ge 
spacecraft material including the effects of hydrogen-helium interactions, inter- 


actions with material defects, defect production by irradiation, and interaction 
of near-surface gas atoms with the materi y to design 
materia 


surface will be neces: 


s resistant to bubble formation, cracking, and blistering. 

Molecular dynamics and den: mulations are appro- 
priate tools to study these effects in silicon, but in order to include realistic 
fluences at high energy in atomistic simulations, one must have access to signif- 
icant computing resources, and artificially enhance the flux by many orders of 
magnitude to reach appropriate fluences in computable model times. Contin- 


uum modeling of bubble formation and bubble dynamics [Blondel et al.(2018) 


offers another approach, but relies on data from extensive HPC-scale MD calcu- 


; functional theory 
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lations. Barring terrestrial experimentation, a simple analysis such as presented 
in this work, however, can be used to select promising shield materials with a 
significant safety margin to compensate for the lack of a complete multiscale 
model. 


4 Discussion 


4.1 Implications on Mission Success 


tions of the previously described phenomena are mission-threatening. 
year journey without an adequate shielding system, hydrogen bub- 
ble formation and subsequent. bursting may induce damage to the spacecraft's 
leading edge which, as the damage continues to make its way into the space 
craft body, will eventually damage components and electronics essential to the 
spacecraft’s performance. Other spacecraft geometries, such as the nee 
spacecraft proposed by Lubin [Lubin(2016)], may help mitigate th 
minimizing the front-facing surface area. If blistering and exfoliation does not 
occur, migration of gas atoms and bubbles into spacecraft components may lead 
to hardware failures, Additionally, over-pressurized bubbles bursting at the sur- 
face will both drastically alter the surface morphology and induce torques on the 
spacecraft, both of which will likely alter its long-term trajectory and orienta- 
tion. Torques induced by bubble rupture must be compensated by the attitude 
em (ACS), which might pose signi ints on the design of 
the ACS and its authority over such perturbations. Moreover, secondary parti- 
cle production and heavy species impacts have the potential to induce damage 
at a greater rate, depth, and severity than the common proton impact. 


4.2 Mitigation Strategies 


Various shielding schemes are proposed as mitigation strategies. Since the space- 
craft is designed to fly edge-on into the ISM, a circumferential shield is proposed 
to protect the ISM-facing edge. The shield must satisfy three essential requisites: 
1) block incoming proton radiation without structural failure due to hydrogen 
bubble formation, 2) stop any incoming heavier species, and 3) mitigate the 
production of secondary particles and their propagation through the spacecraft. 
Since the majority of the expected particle impacts will result in hydrogen build 
up at a localized depth below the surface, the Bragg peak, an effective shield 
would be at least two times as thick as the Bragg peak is deep. 

A favored shielding scheme consists of a sintered material at the H/He Bragg 
peak. The sintered nature of the material makes it granular with a chara 
teristic grain and hole size that enable a high resistance to bubble formation 
[Rossing et al.(1977)]. In a similar scheme, a porous ceramic material is used 
to surround the Bragg peak. The porous nature of the material acts similarly 
to the sintered material in that it makes it resistant to bubble formation and 
allows storage of the implanted hydrogen in its pores. A favorable material for 
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Figure 11: Schematic of a wafer-like spacecraft geometry and four proposed 
shielding schemes. a) Sintered granular material, b) porous ceramic, c) 
self-annealing liquid layer with heat exchanger (yellow), d) metallic powder. 
Other shielding schemes involve plastics, oils, and/or circulating liquids. 


the shield material is carbon due to its higher energy onset of secondary neutron 
produetion than most materials. 

Schemes involving metal powder layers surrounding the Bragg peak are pre- 
dicted to perform in a similar fashion to the sintered material shields due to 
the granular nature of the material. Other shielding schemes involve layers of 
self-annealing metals that can be induced from solid to liquid states and back 
again repeatedly, such as gallium or mercury. In these schemes, incoming pro- 
tons are collected at the Bragg peak while the metal is in a solid state. Then, 
once the rotation of the spacecraft carries a particular circumferential region 
out of the direction of the incoming radiation, the metal is heated until it liq- 
uefies, allowing the implanted hydrogen, in the form of gas atoms or hydrides, 
to diffuse to the edges and escape the spacecraft. Similar schemes involve the 
use of long-chained polymers, plastics, or circulating oils or liquids to mitigate 
bubble formation and store implanted hydrogen. 

More outlandish shielding schemes involve the deflection of incoming charged 
particles via electromagnetic fields about the circumference of the spacecraft. 
However, inherent issues such as high power consumption and strong field am- 
plitudes limit the viability of this approach on a small spacecraft. 


4.3 "Terrestrial Experimental Design 


Terrestrial experimental investigation of the effect of ISM impacts on relativistic 
spacecraft will rely on an appropriate accelerator facility. Such a facility would 
need to meet the following requirements to be informative for this study: 
1. High flux such that ~20 years worth of ISM fluence can be simulated in 
hours, days, or weeks. Since we expect a total fluence of approximately 
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Figure 12: Proposed experimental setup for radiation studies using a linear 
accelerator or cyclotron. Top: The target is to be fabricated in a cylindrical 
shape so as to study implantation over a continuous range of incident angles. 

Bottom: In the ideal scheme, the target mount would house multiple targets of 

various materials and the accelerator facility would have the capability to 
control the temperature of the target. 


4x10? m-?, we would require a flux on the order of 109 m-?s-! to 
achieve the target fluence with an exposure time of roughly 1.5 months. 


2. Ability to defocus the beam such that the ion flux can span the forward- 
facing cross-section of the target, investigating the possibility of blistering 
occurring on surfaces perpendicular to the direction of exposure. 


e» 


. Cooling capabilities to control target temperature. At laboratory-relevant 
fluxes, we expect significantly higher target temperatures than for a rela- 
tivistic spacecraft. 


4. Ability to implant hydrogen and helium at constant velocity. Simultaneous 
exposure may not be possible, so careful alternating exposures may be 
required. 


There do exist facilities that meet many of these requirements. Material di- 
agnostics would include TEM/SEM/AFM of front-facing and perpendicular sur- 
faces to determine surface morphology and characterize any swelling, blistering, 
or exfoliation, and FIB to measure depth-dependent composition to determine 
diffusion coefficients. 

Nuclear activation of the target due to proton irradiation is an additional 
concern for the spacecraft, but is also a concern for access to samples soon after 
experimental irradiations. A future paper will quantitatively address this issue. 
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5 Conclusion 


In this work, we present the potentially mission critical phenomenon of damage 
to relativistic spacecraft by gas accumulation. Composed primarily of hydrogen 
and helium, exposure to the ISM at relativistic velocities will result in gas atoms 
accumulating many atomic layers below the front-facing surfaces of the space- 
craft. Accumulation of hydrogen and helium in solids is associated with pro- 
ach as bubble formation, swelling, blistering, and exfoliation. Since the 
exposure will span the entire forward edge of the spacecraft, blistering may occur 
on surfaces perpendicular to the direction of travel. Using a zero-dimensional 
model of blistering onset, we have performed preliminary calculations in general 
and for copper, an example metallic shield material. For shield materials with 
a yield strength well in excess of 10 GPa, assuming an ISM number density of 
1.0 cm-?, this model predicts gas concentrations will remain below the critical 
concentration for blistering. Using an estimate of the appropriate yield strength 


ces; 


ing synergistic effects of mixed hydrogen and he 
this zero-dimensional model suggests that copper, an example metallic 
material, will blister from hydrogen implantation only if the ISM number density 


as hydrogen-helium mixing and the proximity of stopped gas ions to surfaces 
perpendicular to the direction of exposure may lower the critical concentration 
enough for blistering to occur beyond the capabilities of the zero-dimensional 
model to determine. 

Regardless of whether blistering occurs, the relatively high fluence suggests 
bubbles will form in the spacecraft shield. While no single theoretical framework 
exists to quantitatively explore bubble formation onset and the effe ; 
bubble formation in fusion PFCs provides a qualitative understanding of this 
phenomenon. Bubble formation leads to swelling and a decrease in the struc- 
tural integrity of the spacecraft shield. Experiments or atomistic simulations 
of gas implantation and bubble formation with a wide-beam geometry need to 
be performed to make informed engineering decisions on shield materials for 
relativistic spacecraft. 

Our zero-dimensional analysis suggests a number of strategies to mitigate 
bubble formation and blistering. First, increasing the diffusion coefficient, either 
through heating or material choice, will decrease the local gas concentration, 
preventing bubble formation and blistering. Second, a material that can accom- 
modate large amounts of trapped gas, such as granular and sintered materials, 
will be resistant to blistering and swelling. Third, amorphous materials or liq- 
uid metals are particularly resistant to the effects of gas accumulation. Fourth 
geometric optimization of the spacecraft such as high aspect ratio “needle-lik 
designs minimize the effective exposure area. Without accounting for the effects 
of bubble formation and potential blistering, spacecraft shields may erode or fail 
significantly sooner than instantaneous damage calculations would suggest. 

Total particle fluence depends only on the travel distance and ISM number 


s there 
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y. There exists no feasible mechanism, for low mas: 
ISM particl speeds, so any 
must have a mitigation strategy to deal with gas accumulation. Determining 
the ISM number density along the path of travel for an interstellar mission will 
be crucial for determining the fluence to which a relativis ipacecraft will be 
exposed. Blistering and exfoliation will erode the surface of any future 1 
to destinations significantly further than a Centauri without appropriate miti- 
gating strategies. Regions of particularly high ISM number density will prove 
especially challenging to shield materials. Supercritical fluences for solids will 
require advanced mitigation strategies, such as self-healing, sintered, granular, 
or liquid metal shields that can accommodate the relentless accumulation of 
ISM particles incident upon future interstellar spacecraft traveling beyond a 
Centauri. 


spacecraft, by which 
interstellar mi 


could be deflected at relativi: 


ions 


6 Diffusion 


Effective diffusion coefficients in materials with a significant number of trapping 
sites can be estimated using an equilibrium trapping model [Oriani(1970)]. For a 
single kind of trapping site, the effective diffusion coefficient Deg is determined 
from the lattice diffusion coefficient, Dr, = Do exp(—Es/KkT), the number of 
trapping sites per lattice site, cr/cz, the binding energy of the trapping site Ey, 
and the temperature T: 


Do exp (—Ea/kT) (6) 
1+ (er/er) exp (Es/kT) 

To demonstrate the effect of radiation-induced vacancies on the effective dif- 
fusion coefficient for implanted gas atoms in the shield of an interstellar space- 
craft, an example calculation for hydrogen incident upon a tungsten shield of 
a spacecraft travelling at 0.2c is shown. Tungsten was chosen because it has a 
relatively high hydrogen mobility, and 0.2c was chosen because it has the low- 
est vacancy density production per fluence rate, R = 4 x 1075 vacancies per 
Angstrom-ion, of the three speeds considered in this work. The binding energy, 
Ey, for hydrogen trapping in tungsten vacancies is 1.4 eV [?]. E; and Do for 
hydrogen diffusion in pure tungsten are 0.21 eV and 5 x 1075 m?/s respectively 
[2]. The straggle of hydrogen in tungsten at 0.2c is 20.72 jim, from Table ??. 
Rearranging Equation [6] into Equation [7] the fluence for which enough vacan- 
cies are produced by the incident hydrogen at the effective diffusion length, 
Leg = v Dert, is negligible (that is, below 1 pm) at 300 K is ® = 3.2 x 10? 
cm7?, a condition reached almost instantaneously when traveling at relativistic 
speeds through the ISM. Based on the negligible diffusion due to the effect of 
radiation-induced vacancies alone, diffusion will be considered negligible for the 
energies relevant to this study. 


Des = 


(Do/L2q) exp (—Ea/kT)t — 1 


z= Rexp (Ey /kT) 


(7) 
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